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ABSTRACT 


Methods of generating and properties of Walsh functions 
are discussed with emphasis on discrete Walsh functions. 
Hardware and software versions of two methods of Walsh func- 
tion generation are shown. The Walsh transform is explored 
and hardware and software realizations of a fast Walsh trans- 
form algorithm, particularly applicable to hardware, is in- 
troduced. Computer programs containing the various algo- 
mithms are ineluded. A disewssion of the feasibility of the 
maemo. Walsh fumetions to produce a digital filter with a 


desired frequency response is presented. 
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IL. HISTORICAL BACKGROUND AND SCOPE OF Sine > 


A pair of functions are Said to be orthogeonat= 71 thee. 
tesral of the product of those” functions Over alt opace 
vanishes. Systems of orthogonal functions have been avail- 
able £Oe Many years. At the beginning of the twentieth cen- 
tury, there were many well known and useful orthogonal 
MmMmrecklons: the trigonometrile= functions which occur ime ourier 
series; orthogonal polynomials such as those of Legendre, 
Hermite, and LaGuerre; Bessel's functions; the Sturm-Lioville 
Series; and other special functions. The general theory in- 
emudaing all such eee of functions was essentially devel- 
oped prior to 1910 and it was shown that a set of such 
functions is by no means exceptional but can be produced at 
peel 

Alfred Harr [Ref. 10] formulated an orthogonal system of 
mrmetions in 1909 (now called the Harr funetions), which 
essentially take on only two values. An exemplary set of 
Harr functions is shown in Figure 1. A formal expansion of 
mimearbitrary continuous function™in “terms “of Harr = tunections 
Sonverges uniformiy to the given function. This combination 
of convergence and bi-~-valued properties had not been pos- 
sessed by any other orthogonal set known to that time. 

Independently, and at about the same time in early 1922, 
H. Rademacher [Ref. 20] "and JVE. Walsh (Retr. 24) cache 
bited a set of orthogonal functions taking on the values + 1 


on the interval [0,1]. The Rademacher set, [Figure 2], was 





not complete as other non-trivial fumetions €xisS reece 
orthogonal to all of the members of the set. The Walsh 
functions, [Figure 3], are a complete set and) exhib cocina fan. 
fetes to the trieonometrtemiunctions Sum manyee: thelr peop. 
ties. The Walsh functions are discussed in detail in the . 
remainder of this thesis. 

Very little work was done with these various sets of bi- 
wabuwed LUncEPons until the digital computer became a useful 
tool of the scientific community. Because of the dual value 
characteristics of the Harr, Rademacher, and Walsh functions, 
the applicability to the binary configured machines seemed 
natural. A substantial amount of effort has been directed 
into application areas during the past several years and pro- 
gress has been made in many fields which in’the past have 
only been associated with trigonometric functions. Some of 
the Walsh function application areas of interest include: 
multiplexing, pattern recognition, image processing, electro- 
meatier ic radiation, and digital filtering. This paper wiit 
primarily discuss the area of digital filtering although many 
Mieeehe cresults are applicable to other areas of interest. 

Previously, analysis of digital filters and implementa- 
Pron has been restrictéd to Laplace and Fourier methods and 
hence to trigonometric relationships. The implementation of 
a digital filter based on Walsh functions offers definite 
advantages unique to a two valued function but also has dis- 
tinct disadvantages inherent in the use of a set of discon- 


Einuous functions to describe a continuows | tune eon ELTLToOrts 





expended thus far in the area of Walsh function digital 
filtering appear»to wbesprimarily comcernedgyo ti scbtara aa. 
fhe expansion of a given funcEion in terms vor Seneavaran 
functions and the recombination of the weighted Walsh func- 
m@eons to obtain the filtered and reconstituted funcrton sein 
his thesis the manipulation of the Walsh function represen= 
ation of the given function will be treated and compared £Eo 
Em@andard Fourier methods. Dt will be shown that in general 
there is no advantage to wsang Walsh functors for jdi¢a cal 
filtering when a specific frequency response is desired, but 
that for certain frequency responses, Walsh functions do 
Zive a decided advantage over the Fourier methods. 

Section II presents some of the properties of Walsh func- 
mons including a modified notation method which, it is 
peuwmsevyed, is clearer and easier to usé in mathematical re= 
iationships. 

sections III and IV show methods of Walsh function gen- 
eration and Walsh transform techniques respectively. Algo- 
rithms are presented which are adaptable to either hardware 
or software and will either generate a set of Walsh functions 
Or will take the Walsh transform of a given function. Serial 
maeparallel input of a sampled function is provided for. 

section VY considers the wse Of Walch functions seo mpmocdmac 
a digital filter with a desired frequency response. The fre- 
quency spectrum of discrete Walsh fumetions) 15 ert 7edewang 


it is shown that in general it is as complex to combine these 





to achieye a desired frequency spectrum@output nage teem 


case in the Fourier domain. 





IL. PROPERTIES .OF WALSH UNC REONs 


Walsh functions can be defined oyer the interval [0,1], 


ieene Lollowine mMapmer ; 
ci, O04t<l1 


DRC 0, otherwise 
Kan 
W(2k+p,t) = WCk,2t)+(-1) W (kee) (1I-1) 
pea 0,4 
k= 0,1, 22 


When Walsh functions are used for digital filtering, the 
input function is normally sampled at some rate as are the 
Memesh functions. En order £o discriminate barwean the dis- 
crete (sampled) and the continuous functions and for ease in 
Mathematical manipulation, the following notation is used? 


Continiuowsc CS ted NOt con. 


W(k,t) 
Tee Ob 2 ane index Of Walsh fumenrton 
Of£tel1 independent variable 


Dascrete funckzon notation: 


Ko 0.124.246 N=) tide womeWwash funce2 om 
N = a power of 2 Gasdeyr of Set of Wallsh@inmetrons 
j = 0,1,2,..., N-1l independent variable 


The discrete Walsh fumetions are formed trom Pnemeome: 1— 
uous functions by sampling Wtk,t) at t = Gunn Ciese 4. fem 
the right). Figure 4 shows the set of order 16 discrete 


Walsh functions arranged as a 16 x leumarri 
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Several other forms of notation have beeneucede ofan 2 
functions. Some of the most common of these are tabulated in 


Table 1 for comparison) andm@ercess Trerenencce pimp ascen 


Discrete Continuous Walsh Corrington  Harmuth 
aaa os Pere aye ~ 54045 
[Ref. 9] 
We (0,35) W(O,t) >t) Vogo 6) Cal (0,9) 
We (1,4) W(1,t) >, (t) Yoo bt? Sal(1, 8) 
We (204) w(2,t) @ Gt) vy Ce) ead (2, @) 
We (3.4) wi3,t) F(t) Wy g(t) 8a (2,8) 
We (454) wast) 6 Ce) yg le) Cad (2, 8) 
We (555) wOS,t) 629 (tE) Wy Ce) 8a (3, 8) 
We (654) wi6,t) 6 BPCe) yg, CE) a (3,8) 
We (743) w7,t) S(t) Wyo lt) Sa (4,8) 
Table 1. Notations @#on, Walsh Fumctions 


These functions form a complete, orthonormal set over 


Meeminterval [0,1). The impltecation of this property vs Celiat 
any absolutely integrable function defined over the interval 


can be expanded in a series of Walsh functions analogous to 


mem OuUridGr E€xpansion of such a function, 


ier 1s: N-1 
E(t) = 2 a, Wk, t) (II-2) 
k=0 
where: 
1 
ay = Jf £Ct) Wk,t) dt (II-3) 
0 


Assuming f£(jt) is a sampled representation of f(t), 


then a, = Q for 


where the sample period, Tt = 1/N; . 
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k IN aide Gina teens (ITL=2) and €1i-3)  vwedwecwec. 


N-1 
fi) oe ay Wy Ces 5) (LI-4) 
j=0 
and 
N-1 
ay may NN Ae EC) eon (II-5) 
j=0 : 


Lf the function £ C(t) soree () 7) otsedetinedsover (rnc int ciae 
mat [0.17], a simple change if scale, t, = t/t willealilow Sthe 
use of the above equations to find the Walsh representation 
oii Ae fimecetonyoever [0 bi- 

Sequency is defined as one-half the number of zero cross- 
ings per unit time. The number of Zero crossings and hence 


the sequency of a given Walsh function can be obtained by: 


_ ee (nt? 


5 (II-6) 


where zis the sequency 


and k is the index of the 
Waish function 


The set,ia,}, therefore represents the sequency spectrum 
mee. t) in the same sense that a Sét of Fourier cockElectents 
represents a frequency spectrum. It should be noted that 
the Walsh functions defined in this thesis are ordered ac- 
SBerding to increasing séquency. This is the definition as 
used by Walsh in his original paper [Ref.24] and is equiva- 
dent to that used by Harmuth [Ref. 9]. Another method of 
defining Walsh functions is as products of Rademacher func- 
Eions. When defined in this way, they are mot ordered 
according to increasing sequency but are what is termed 


binary ordered. Appendix A shows this method of defining 


ty 





Walsh functions and derives algorithms £07 §eenvet Pigeon 
one method of ordering to the otheu. Sine (on eea pale, 
ordering of Walsh seems, to be best suited tothe spi mpgce ca. 
communications engineering. 

Another property of Walsh functions can be seen from 


their internal symmetries. A product of two Walsh functions, 


if they are both even or both odd about 7 will be even. If 
Ome is odd and one is even, their product will be odd. This 


is also true with the symmetries about a, -, € tC... ot nema ds 
product function is completely described by its symmetry 
Pueperties about the points Same 1 1 a he ee atest See 
concluded that this product is also a Walsh function. More- 
over, the index of the product is the modulo 2 sum of the 
mmdices of the two factors. Thus: 
WCk,t)*Wm,t) = WCk@n,t) 
where @® indicates addition modulo 2. 


Ie ren C245) Wo be W(10,,t)-W(11,,t) 


W(10,@11,,¢t) 


I 


WC(O1,,t) 


WG we) 


I 


muce the reciprocal of a Walsh function is the same function, 
WC te Wii e) 
the Walsh functions form a closed set under multiplication 
Sr diyision-. 
aA usetulyproper by #0 f mds sicmoeuse Walsh functions is that 
Wy (Gok) = Wy Ck, 4), 


that is the matrix representation, Figure 4,15 symmeecmce 


oe. 





A Hademard matrix is an N x N matrix, each of whose ele- 
ments are tl or -1 and Such that Ehe Vews wo bet ne emcee 
orthogonal to each other using the ordinary vector inner pro- 
duct. Since the discrete Walsh function representation con- 


Bemms tQq this definitiomey it is aesymmetric wim amarde mae, 


If H denotes the matrix, then: 


H*H = NI where Ale iseH Eranspose 
N is the order and 


lis the identity @atrix 
but H is symmetric, so H = hee and therefore, HH = NI. This 


important property will be used in Section IV, The Walsh 


meanstorm. 
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Til. GENERATION OF WALSH FUNG ro 


Walsh functions can be generated in seyeral different 
ways, all of which depend upon the mathemauical=propemeres at 
the functions as discussed in Section II. Figure 5 shows a 
hardware realization of a Walsh function generator designed 
by Harmuth [Ref. 9]. This generates Walsh functions as pro- 
ducts of Rademacher functions as discussed in Appendix A. 

Swick [Ref.22 and 5] showed that by using the properties 
of symmetry, a particular Walsh function can be written down 
merectly from deft t0 Fight without making use GCE any preced— 
ing function. This method is to expand the index of the 
desired function into binary. 

eer Kk = Dee Db wie N= 2 


n-l nee 
and let W(k,t) = 1 in the interval [0,1/N]. Then for j =n, 


> 


=) 


oll, Sop Ue Ty Seen aya iCk, 2) ts edlilicvan) ahowe 2 2 
when it is considered as a function over the interval 
oe Figure 6 shows the generation of W(11,t) by 
Mies method. Eleven in binary is L011 and nm = 4. Weteing 
feelte.t) = 1 in the interval [0,1/16], then b, = 1 indicates 
ameodd function about 1/16” on the interyal [0,1/3]- b. = 0 
midiecates an even function about 1/8 over the interval 

mo, 1/4): b, = 1 indicates am odd functwon about |) 4ecver 

MO. /2)- and by = l indicates an “odd Funct von sey ct ge cme meee 
moreerval” (Oe). 


Peterson, [Ref.18]}, devised a method Similar towcwec come 


Which the index is expanded into reflected binary (Grey code). 


digs ‘ 





Hei Kalk ++ 81> W(k,t) = 1 in Uehtewinter, aioe 


Sy 5y-1° 
and then for j = m, n-1,:.. Wk, teehee ore 2 

ios po le is equalito (25) en - yee in the interval 

(i, 2a Referming again to Figure 6 f@aamiehe generation wo 
molt), and using the algorithm from Appemawers = toeecon ven 
miom from binary Go reflected binary, ll an) xetlected man sn 
meelll0. Letting WCll,t) ="I@an the interval” [0577 16) eeemem 
(-1)®4 = -1 and W(1l1,t) in the interval [1/16,1/8] equals 
Mmacerlt) an the @mterval [OS1/lo).. Similarly, woes aa. 
mwee.)/4)] equals =W(11,t) in [0,1/8)2 WOll,f), 1/4 tay 
eemals -W(1l1,t) for 0 t 1/2; and the function from © = W772 
mom = 1 is equal to WC1l1,£f) over [0,1/2]. 

Lt should be noted that by using Swick's method but ex- 
panding the index in reflected binary or by: using Peterson's 
method with the index expanded in binary, the Walsh functions 
generated are in binary order instead of sequency order. 

Figures 7 and 8 show hardware versions of a fast Walsh 
transform algorithm which will be discussed in further detail 
mieoection IV. A single pulse input into the device of 
Figure 7 or a pulse shifted through the input register of 
Pigure 8, will produce as an output the Walsh functions in- 
stead of the sequency spectral components. 

Computer program 1 utilizes a set of algorithms [ketsli- 
12, and 13] to form a Hadamard matrix of Walsh functions in 
either binary or sequency order. These algorithms compute 
the value of the Walsh function at each point without, verter. 


ence to previous points or previous functions. slo. tome 
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relatively slow method of function generation since it does 
not take advantage of the symmetry properties cf (thes var-n 
functions . 

Computer program2 also generates a Hadamard matrix in 
sequency order using an algorithm based on Peterson's method 


of Walsh function generation. 


17 





IY. THE WALSH TRANSFORM 


It was shown in Secitomelimeenat a functions canmspenc.. 
pressed in a series based on Walsh functions analogous to a 


bOourier serieS. olin this series, a set of coefficients, 


{a 


Ko? represented the sequency spectrum of the given function. 


The discrete or sampled form of the Walsh transpose became: 
N-1 


a. = 1/N © £(47)°W. (k,j) Cy) 
k 4=0 N 


Menthe samples of £(j1T) are arranged in a column matrix, 
oO) 


fey) 
£O7) 


x 
fl 


£{ (N-1)T] 
and premultiplied by the Hadamard representation of the Walsh 
mimetions, equation (I1V-1) is realized in matrix form as: 

W = 1/N HF | (IV-2) 
waere W is a column matrix of the spectral coefficients. 
Because of this property, the finite Walsh transform is often 
called the Hadamard transform. 

From the properties of Has described ine section. Ui 
ean be readily seen that: 
HW = 1/N HHF = NI/N F = F (IV-3) 
or that with the exception of a SORE One multiplier, 1/N, the 
Same procedure is used to obtain both the Erans orm deere 


inverse transform. 
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Some of the properties of the finite Walch ewan 
Ref. 15), aresas) Gottiows: 
if fC ee, £ 
and g (jt) (g) 
ee erans hommepaiirs, then: 
(1.) The transform is Peinear 
ai (j tee be e (11) 7a CE) th We) 


mere real constants a and b. 


(2) The coefficient of index 0 is the sample mean of f(j) 


N-1 
ape = iy7N 2 £ Cit) 
j20 
Mmoeical or dyadic convolution of two sequemees £()) and 2()) 
is defined as: 
_ Nel 
mC jT)SML/N OL OF CK Ue CK) 7} (IV-4) 
k=0 
where @® indicates modulo’ 2 adda Cajon 
(Sy 
AC) = Gielen ed) eee GE mca (Sy) 
where g indicates logical convolution 
2 
maa (4) £(j7T) O£ (JT) +——-WO) 
This is the analog of the Wiener-Khinchin theorem of the 
Mmeequency domain. Ene Walsh power Spectrum 1S ene faneeee 
feish transform of the Poercal autocerrel a: Homme mee 
Computer” program. 2. We wdutzes the Hadamard matrix multi- 


plication method to give both the transform or inyerse 


transform of an input. This method is stralgenmt fPorwavdeou 
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time requirements can be decreased considerably by using fast 
transform algorithms much Jike the fast Fourter trans ou 

Several approaches exist for fast Walsh transform imple- 
mentation. One method described by Shanks [Ref. 21] parallels 
the derivation by Cooley and Tukey for the fast Fourier trans- 
form, but even faster methods are available [Ref. 2, 23, 
and 25] which make use of the internal symmetry properties 
rmethe Walsh functions. 

A useful method is closely associated with Peterson's 
method of Walsh function generation. Instead of the reflected 
binary expansion of the index indicating whether the inter- 
val is repeated with or without a sign change, the index 
mide cates whether the interval is added or subtracted with 
Eee adjacent interval. For example, for the second spectral 


component of an eighth order system: 


fof totazt, fiefs 


(reflected binary) 


Ee) 
eno tl 

8; = 9, €& = 1, 8, = 1 
and a, = 1/8(£ t+£,)-(£,+f,) - [(f£,t+£,)-(£,+f,)] 

Prd ft ff 4 

on. ae ae ee 
Higure / depicts a hardware realization of this type system 
hoy an eighth order spectuum,. 

Using this or other fast Walsh transform algorithms 

reduces the number of operations required to NlogoN real 


additions (and/or subtractions) as compared to Nlog,N complex 
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multiplications and Nlog.N complex additions as required by 
the fast Fourier transtormeatcorittime- It is therefore much 
faster to go from the time to the sequence, @domai nel ao 
the time to the frequency domain. Also the inverse Walsh 
transform can be performed with the same hardware or software 
as the transform and is therefore much more economical if the 
cost is calculated in time, computer memory or hardware. 

Computer program 3 and Figure 8 depict software and hard- 
ware realization of another fast Walsh transform method. In 
this method, sums and differences of the sample values are 
formed according to a pattern. The sequency spectrum coeffi- 
cients are the output and in the hardware version, the se- 
quency corresponding to the coefficient can be found by trac- 
ime the most direct meat from anyone of the sample input 
Peimts to the desired output. When an addition is passed 
through,a zero is written down; when a subtraction is passed 
through, a one is written down starting at the higher and pro- 
ceeding to the lower digits. When the output terminal is 
meached, the digits written down are the refiected binary 
mepresentation of the index of the coefficient. 

The software version, represented by subroutine FASWAL 
of computer program 3 requires the subroutine eee CONVRT 
mo convert the coefiiecients to thessequency crder.,  sliece 
coefficients are in a reversed reflected binary order much 
as the discrete fast Fourier transtorm coro t put cera cme 


GBeversSedwpinary soader , emer oO |. 


Zak 





To find the corresponding indicies as CONYRT does where 
W(m) = Tk) 
W(m) is the Walsh coefficient in sequency order 


T(k) is the coefficient in reversed reflected 
Dina vy sio1 diem 


(1) expand m into binary 


(2) convert to reflected binary 


m = one Seay 


(3) reverse order of digits 


k = G185--°8 


n 


(4) use this as the binary representation of k 


Figure / shows another version of this same fast trans-— 
merm algorithm. In this case, the input is applied serially 
and after N input samples, the output (in parallel) will be 
the Walsh transform or the inverse transform if the input is 
me, Sequency ordered spectral coefficients. Note that as 
mor Fieure 8, the output is in reversed reflected binary 
Sraer. 

Either of these realizations have advantages when hard- 
mane implemented since the Ordering ofp mote pun scalenc 
predetermined and hard-wired so the reordering is not neces- 
Sary each time the system 1s used.  Noittewalsowenat ec meter 


these units can be used to generate Walsh functions. The one 


Ze 





in Figure 7 by putting a unit pulse sin tose iew inp ie sande 
one in Figure 8 by shitting a unat pulse down tihwowe ee 


impuG register. 
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V.. DIGITAL FILTERS CONS PDE Reyes 


One of the traditional methods of digital £ileeriiceen 
the frequency domain is by the use of the convolution theorem. 
Mnat is, an input function is transformed into, tne ftrequence, 
qomain usinp either the Fourier, Laplace, ox "Z” tranistiorm. 
thas is multiplied by a weighting function, the transform of 
the impulse response of a desired filter, and the inverse 
transform is taken, yielding the desired output. The logical 
convolution theorem does not produce this result. 

Much has been written about the transition to and from 
the Walsh or sequency domain for the purpose of digital fil- 
tering but very little material can be found concerning what 
can be done while in the sequency domain to obtain the desired 
response from Bey Aen cer in the time domain. 

It is possible to construct a filter with any desired 
sequency response. There is a rough correlation between se- 
quency and frequency responses of filters but only to the 
point that the higher sequencies contain more of the higher 
frequencies than the lower sequencies do. For example, low 
pass sequency filtering at a sequency cutoff of one-half the 
bughest séequency, the sampling sequency, is equivalent te 
halving the sampling rate of the input function. 

To take advantage of the high speed and other economical 
properties of the Walsh functions, it would be desirable to 
be able to obtain a desired frequency response of a digital 


filter which operates in the sequency domain. 
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The first step in the investigation of Elis posotp 1a 
will be to find the £requency spectrum of the discrete en 
functions. (The frequency spectrum of the continuous Walsh 
functions is shown in Appendix B.) 

The Z-transform of a Walsh function can be written as; 


N-1 


ZIW (k,i)} = 2 Wy(k,j) 2 (V-1) 
rod 
moting that for N = 1, 
Zit OmO) ae = a (v-2) 


-1 
Z{W, (0,3) ] = Mapes 
al 
Z2(W,(1,3)] eae C¥e2) 
the following recursive relationship is postulated: 
; , i+ N/2 
ZIW, (2k+p,4)] = ZEW, 75 (k,5) )fat¢-2)3*P 28/21 (va) 
N N/2 
oe Oe Ne ee es ee 
Wee PO ale? a ete Noe N = 2,4,;8, 


ints can be induced by the following argument for N a power 
of 2, say N = aaa Expresisane Kk of equation (V—4) in biimary 
elie) S 


os iter Saad (V-5) 


Bue binary representation of k can “be converted to ret Peerca 
binary (see Appendix A) as; 
iu Ce etna al 


® babe eee 
(V-6) 


Zo 





Noting that, 


[2-k], = b Dis seamis b,0 (V-7) 
and lie = Bi 8y-1 61°81 °y (V-8) 
then, BEM a yoo ola Tl (V-9) 
and 
mele = g 85-1757 8z> (V-10) 
where rep 50% ee ae 
QO for bs = 1 


From Peterson's method of Walsh function generation as 
Gescussed in Section IIL, the higher order function is com=- 
posed of the lower order function over the interval [0,1/2] 
and it or (-1)°1, (for 25} and 2-341, respectiveryore 
memes the lower order functiome over the interval [1/2.1)% 
Since ane = (ied and ys = (-1)°1, equation (V-4) is 
established. 

The Z-transform of a Walsh function will therefore be of 
Ene form, 


eS Cem 


ZUWy(k,3)] =(1t277) (it27°) te")... (tz 
As can be seen from equations (V-4,7,8,9 and 10), the 
sequence of signs is determined by the reflected binary code 


Seethe index, k; that is, 


: _ SG eo . 
atWy (kG) = TP (1ec-1) "2°73 (V-12) 
i eer 

r= ®[1og, (N/n)-1] where g is the subscripted digit of 


the reflected binary representation 


Of. ke 
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evaluating, 





ee ee eo es ee ST 
z=e 
1 = sample (pemiog 
W = radian frequency 
= LmWT -im)T -imWT 
= 2 (ep e wear eae 
Z 
= 2Cos(mwt)e ~p— (V-13) 
4 
letting x = = — ne i = 1/t = sampling frequency 
Ss 
Bnien, 
142ml 2 . = 2Cos(xm)e ~~~ (V-14) 
si {ub} 96 —— Z 
zZ=e 2 
paepilarly, 
-m/2 a =sel pu T S 
l-z eeivt fer any 5 + 5 (V-15) 





Therefore the frequency spectrum of a discrete Walsh function 
fimcaampliy a product of a group of sine and cosine functions 
with a linear phase response and can be written down directly 
ion expansion of the index, k, Me reflected binary. For 
Wy Oko 4d); expanding k into reflected binary eaves, 

k Be .s40° oa = log, (N) 
then the magnitude of the frequency aneenat component is 
Bayen as, 


Pees) || imei sO COs era Gon we 


eee Costeemetor 8 
P Shalgel '@ Mae oye ae 


x) | CVG 


I 


and the phase is eilvenepy, 


/¥lilty ke 40} = ~(N-1)xtm7/2 


Za 





where 


Bor example, thestreqwency sopecterumsor Wa, (14,4) would be, 


10 2 
N 


k = 14 = 01110, = Be. 


II 
wo 
Nh 
| 

I} 
in 


8. 8),858584 = 01001 


JFIW,, (14,4)]| 2? | cos (x) Sin(2x) Cos (4x) Cos (8x) Sin(16x) | 


and the phase is, 
Z FIW,,(14,35)] 


Figure 10 shows the amplitude of the frequency spectrum of 


= Bix rece = Sar TW 


oe 
| 
| 


@ne first eight Walsh functions of order 16 plotted on the 
interval [O#f/f <.5]. Note that 0.5 f/f. is the maximum or 
Nyquist frequency for digital filter applications. eri 
general it is aa Temeeciary to plot the first N/2 spectrums, 
since by symmetry, the N/2 remaining spectrums, (Figure 11), 
can be obtained by reversing the f/f | axis of the lower order 
frequency spectrums. 

One method often used for digital filtering invthe fre- 
micncy dOmMain is to use digital resonators in the form o£ 
eomp filters with zeros introduced to give the desired spec= 
trum response. A Walsh function resonator can be defined as 
omctgital “cireult whose pulse response 25 cag some ure como 
me Walsh resonator can be constructed as showman fie cee 
Rere, the sign on each input to the N adnput “Summerets gemen- 
Mined by the sign of the particular Walsiesumerton 92 eae 


POLNEs Ju. ek Wa, (14,3) resonator is shown in Figure 13. 
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The z-=transtorun of W.Ck,j), equation (Y-11), leads to 
another and simpler realization of the Walsh resonator. thee 
method, shown in Figure 14, uses cascaded sections of lower 
order filters. The Wao 4,4) pons ioe of this type is 
shown in Figure 15. This is the method which led to the 
Peeelopment of the fast Walsh transform algorithm discussed 
in Section IV. The hardware realization of that algorithm 
shown in Figure 7 is simply a full set of N Walsh resonators. 

Referring again to Figures 10 and 11, it can be seen that 
various frequency responses can be obtained by cascading and 
paralleling various combinations of Walsh resonators. There 
is no method known hee will give the classic filter frequency 
responses such as Butterworth, Chebychev, or elliptic using 
Maas method. A function minimization of the mean Square error 
between a desired frequency response and the frequency res- 
ponse of cascaded/parallel Walsh resonators can give the 
necessary combination whose frequency response is:‘as close 
to the desired response as is possible with a given number of 
resonators. This method will lead to results comparable to 
that of using a digital comb filter and resonators in the 
frequency domain, [Ref. 8]. 

As an example, a filter is desired with a narrow band- 
Width at approximately f/f. = .25, using three Walsh resona- 
tors. Using the Wi 672k) resonator as a starting approxima 
tion, it can be seen that the bandwidth can be decreased and 
sidelobe magnitudes decreased by cascading it with notch 


iiiters "cOnmstructed. 1 soul Wig 665k) and Wig 692k) resonators. 


oo 





The notch filters are simply resonators noce soley tte oa. 
tracted from the input in a feed Sion vwardgmanger leo 
Made possible by the Pimear = pmase Chae t cm iet Neo Oimrene. 
Walsh resonators. 

Figure 16 shows a block diagram of such a filter and 
Figure 17 shows the frequency response of the filter. The 
mab bandwidth is about O4f/5 | and the mearesteertde lobe is 
down approximately l16db from the main lobe. Neither of the 
diagrams takes into account the delay errors which would 
probably be introduced by such a system but these errors can 
be eliminated much in the same way as they are in frequency 
Gomdin digitai™filters. Filters with higher Q's and better 
sidelobe suppression can be obtained by using higher order 


Walsh resonators and/or more of them in the system. 
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YI. RESULTS AND CONCLUSIONS 


Many computer runs Of@bethigan oneline andigit —liae snare 
Moing the Walsh transforms for digital filtering Haye been 
made with vatiews types of Mmputsmandsfalremeresponscese Ai ter, 
this work, it was concluded that the use of Walsh functions 
immnot thesabsolume answer tomdigitalgfilterine om gto tlie 
Many other fields of application being investigated. However, 
Eaey do Annies have definite applications to assist in 
overcoming certain problems. 

The areas, other than multiplexing [Ref. 4, 14, 15 and 
17], where greater benefits may be obtained as compared to 
Beeconventtonal methods appear to be where either the input 
Omoutput function is of a discontinuous nature. One area 
where this has been put to use is in the field of image pro- 
cessing [Ref. 1, 3 and 19]. ‘Here the 2npus Lun et ron ener — 
mally two dimensional, may have noise bursts included in it 
which may be attenuated by low pass sequency filtering or may 
have low contrast areas which can be enhanced by high pass 
mayuency filtering. In the first case, the input has discon- 
Poouirties, the noise bursts; in the second, the output re= 
murires the presence of discontinuities, the image edges. TE 
has been shown that digital image processing by multi-dimen- 
Sional Walsh methods gives essentially equivalent results to 
Fourler methods but with the speed increase inherent in the 


sequency domain. 
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Image processing and the field of pattern fecoen ein 


[Ref. 6] have much jin common and again if the pattern to be 


processed is of a discontinuous form, it is believed that 


Walsh functions may again be superior to Fourier methods. 


Another signal which has characteristics that lend 


themselves to Walsh functions is one of the "infinitely" 


clipped type. Infinitely clipped signals also only take on 


two values and can be represented by Walsh functions easily. 


The use of Walsh functions in this area, especially for 


mpeeication to infinitely clipped 


speech processors and 


woceders should be investigated further. 


Walsh functions will probably 
in the fields of image processing 
the exception of very specialized 


Mera lter sinusoids will probably 


gain even preater status 
and multiplexing but with 
applications, their use 


remain quite limited. 





AE ye eee 


Walsh Functions as Products OF @hagemcn aime 1 Oates 


The Rademacher@it#ictions can be@defined as Hollows: 


até 
RO; Vee ee 
O, otherwise 
1, 04145 
le ae al 
2 
0, otherwise 
then 
R(n,t) = R(n-1,2t)+ R(n-1,2t-1) (A-1) 


Tiere Ze. 4. Meteo, 
mach Rademacher function is a replica of the preceding func— 
tion with the sequency doubled. Referring to Figure 2 and 
Figure 3, the first 5 Rademacher functions are W(0,t), W(1,t), 
Paeont), WC7,t), and WClS,t). That is, the Rademacher func— 
tlons are the Walsh functions numbered ele iy = 0 ee eee 
Beminaing a different method of ordering the Walsh functions, 
W(m,t), and expressing m in binary forn, 


m= g 8,-1°°°8281 


mnen : 
nh 
W(m,t),= 4 | R(j,t) 23 Chee 
Ve! 
J Ati 


ed WC100,,t) RiGee) 


b 
W(S,t),= W(LO1,,¢), =" RO3ob) RG ep 


1 


Se 





Noting that W Cee W Ghee ead WOS,t) = W(6,t), the func- 
tions WGm,t), are notein thessame order as W(ket) -0n lomeend 
the index, m, in Wm, t), which corresponds to the index, k, 
in W(k,t), it is necessary to represent W(k,t) as the product 
of Rademacher functions, remembering that R(i,t)= WiC Ee 
From the product rule, W(p,t)‘W(q,t)= W(p@q,t) and express- 


b the desired 


mao K in its binary digits: k = ee i 1 


relationship can be found. 
From (A-2): 


n n 
W(m,t),= | R(j,t)84i = | KO Sime) a 
j=1 jel 


but gJ-t in binary is of the form Jit--++1l and k mustebe ene 
modulo 2 sum of several factors of this type. This can be 


expressed as: 


> oy SA I reads meee a2 n° By eles eee RET 
ny die actis n=-1 digits 


Adding the right side gives, 


ber Oe g(s 8 Boe Cea B-1 © B29) +++ (8,%---@8,) 
Our , 
ae Boy Ee 198° TB gy IE 
Therefore 
by” DS 
g,° Sor cs! ote 9 eal 
Cr 
els ae Y Drei 
and 
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Thus the binary representation of We Seep OST is the 
reflected binary (Grey code) representation of oc 
ioaat is, 


Ca 


-+b,,¢) a Wg 8a 8a t)y (A-3) 
The following algorithms can be used to convert from binary 


moereflected binary or vice versa. 


Dn ee LY re (=p) 


Mis is illustrated by the following diagram. 


DL eal bo? i 
b * ‘s, =, b b 
n+l on n-l 3 2 
Bn Bn -1B n+ 2 Bo = By 
Dee oe ce aoe 2 2S 
A A A “A 7 
Se Le 7 Dy By 


Figure 9 shows a set of discrete Walsh functions genera- 
ted as the product of Rademacher functions and left in 


Ganary oxder. 
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APE EN DUX« 8 


Frequency Spec hr une (Hota he himeans sonny 


of the Continuous Walsh Functions 


The Fourier transform, F(w), of a function is defined as, 


F(w) = f oe a) dt i= ¥-1 


w= 


W radian frequency 


therefore, 


1 : 
aa) = ere Clee ake 


0 


and by the definition of W(k,t), 
wt -iwt jx : -iwt 
Fi, (w= ie Wik 2 fy det (i see W(k,2t-1l)dt 
0 0 


Or 


Fy, (w)= 1/2 ¥, (w/2) (14 (-1) Se 98) 


pemilarly, 
Fo, ym) = 1/2 F, (w/2)[1-(-1) Se 28 


Therefore, 





FY Cw) = = (1-e ) = 2 ‘sinGe 2 
W 
Dy el -iw 
Fi (w) = 1/2[ 4 Sin(w/4)](l-e =e 
w/2 
= a8 e 48 un Gs 
ee 
Fo (w) = $42 -.°2E sin? (w/8) Sin(w/4) 
Fw) = as e 2X sin® (w/8) Cos Gr 4) 


By induction, a general form of FL Ga) can be established as 


follows: 
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expanding k in reflected binary. 


a en enn’ * 81 
Lieve 
n 
m = - 8 
then, 
£, (w)= 2h" ve 28.sin®(w/2"**). 


n-1 
| | EAGd/2d, Cos acta 
j=l 


Bi 








Fig, i. The firse l6enarr tunecrewoncy Xfi) Ct). 
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Walsh Transform as Hadamard Matrix Product 
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